Magnetodielectric effect and phonon properties of compressively strained EuTiOa thin 

films deposited on LSAT 



S. Kamba/'HV. Goian, 1 M. Orlita, 2 D. Nuzhnyy, 1 J.H. Lee, 3 D.G. Schlom, 3 ' 4 K.Z. Rushchanskii, 5 M. Lezaic, 5 T. 
Birol, 6 C.J. Fennie, 6 P. Gemeiner, 7 B. Dkhil, 7 V. Bovtun, 1 M. Kempa, 1 J. Hlinka, 1 and J. Petzelt 1 

1 Institute of Physics ASCR, Na Slovance 2, 182 21 Prague 8, Czech Republic 
Laboratoire National des Champs Magnetiques Intenses, 
CNRS-UJF-UPS-INSA, 25, avenue des Martyrs, 38042 Grenoble, France 
1 Department of Materials Science and Engineering, 
Cornell University, Ithaca, New York, 14853-1501, USA 
J, Kavli Institute at Cornell for Nanoscale Science, Ithaca, New York 14853, USA 
Peter Grunberg Institut, Quanten-Theorie der Materialien, 
Forschungszentrum Julich GmbH, 52425 Juelich and JARA FIT, Germany 
6 School of Applied and Engineering Physics, Cornell University, Ithaca, New York, 14853, USA 
7 Laboratoire Structures, Proprietes et Modelisation des Solides, 
UMR8580 CNRS-Ecole Centrale Paris, 92295 Chdtenay-Malabry Cedex, France 

(Dated: March 2, 2012) 

Compressively strained epitaxial (001) EuTiOa thin films of tetragonal symmetry have been de- 
posited on (001) (LaA103)o.29-(SrAl 1 / 2 Ta 1 / 2 03)o.7i (LSAT) substrates by reactive molecular-beam 
epitaxy. Enhancement of the Neel temperature by 1 K with 0.9% compressive strain was revealed. 
The polar phonons of the films have been investigated as a function of temperature and magnetic field 
by means of infrared reflectance spectroscopy. All three in-plane polarized infrared active phonons 
show strongly stiffened frequencies compared to bulk EuTiOa in accordance with first principles 
calculations. The phonon frequencies exhibit gradual softening on cooling leading to an increase 
in static permittivity. A new polar phonon with frequency near the TOl soft mode was detected 
below 150 K. The new mode coupled with the TOl mode was assigned as the optical phonon from 
the Brillouin zone edge, which is activated in infrared spectra due to an antiferrodistortive phase 
transition and due to simultaneous presence of polar and/or magnetic nanoclusters. In the anti- 
ferromagnetic phase we have observed a remarkable softening of the lowest-frequency polar phonon 
under an applied magnetic field, which qualitatively agrees with first principles calculations. This 
demonstrates the strong spin-phonon coupling in EuTiC>3, which is responsible for the pronounced 
dependence of its static permittivity on magnetic field in the antiferromagnetic phase. 



PACS numbers: 75.80.+q; 78.30.-j; 63.20.-e 

I. INTRODUCTION 

Multiferroic compounds in which magnetic and elec- 
tric orders coexist are intensively studied due to their 
high potential in magnetoelectric devices as well as due 
to their rich and fascinating physics related to the mag- 
netoelectric effect. Unfortunately, it appears that there 
are only few multiferroics in the nature and their abso- 
lute majority exhibits multiferroic properties deeply be- 
low room temperature* 1 Moreover, most multiferroics ex- 
hibit antiferromagnetic (AFM) order and therefore only 
weak magnetoelectric coupling. 1 For this reason, there 
is an intensive search for new materials with ferromag- 
netic and ferroelectric order, where the magnetoelectric 
coupling could be high. 

Fennie and Rab©2 suggested a new route for prepara- 
tion of multiferroics with a strong magnetoelectric cou- 
pling. They proposed to use a biaxial strain in the thin 
films for induction of the ferroelectric and ferromagnetic 
state in materials which are in the bulk form paraelec- 
tric and AFM. Basic condition for such a material is 
the strong spin-phonon coupling. Fennie and Rabe^ pro- 
posed, based on their first principles calculations, to use 
EuTiOa for such purpose. Recently, Lee et al4 actually 



confirmed the theoretical prediction and revealed ferro- 
electric and ferromagnetic order in the tensile strained 
EuTiOa thin films deposited on DyScOa substrates. 

Bulk EuTiOs is an antiferromagnet with G-type AFM 
order below T/v= 5.3 K. 4 Temperature dependence of its 
permittivity e' exhibits quantum paraelectric behavior, 
i.e. e' increases on cooling and saturates at low tem- 
peratures due to quantum fluctuations, which inhibit 
the creation of long-range ferroelectric order. Recently, 
it was theoretically predicted that EuTiOs has antifer- 
rodistorted ground state££ Experimental studies show 
that at room temperature the crystal structure is cu- 
bic Pm3m perovskitej 7 - which transforms to tetragonal 
IA/mcm structure near 280K.— ~— e'(T) exhibits a sharp 
drop below Tn due to a strong spin-phonon coupling^ A 
large interaction of the magnetic moment with the crys- 
tal lattice was also manifested by 7% increase of e' with 
magnetic field, observed at 2KJ^- Linear magnetoelec- 
tric coupling is forbidden by the symmetry of EuTiOs. 
In spite of this, a strong bi-quadratic magnetoelectric 
coupling was recently reported below T^J^ 

Low-temperature infrared (IR) reflectivity studies of 
EuTiOa ceramics reveal three polar phonons typical for 
cubic perovskitesi 13 ' 14 The frequency of the lowest fre- 
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quency (T01) phonon decreases (the mode softens) on 
cooling, which fully explains the increase of static permit- 
tivity on lowering temperature via the Lyddane-Sachs- 
Teller (LST) relation. Below 100 K the soft TOl mode 
frequency deviates from the classical Cochran law and 
finally levels off below 30 K.— Such temperature depen- 
dence follows the Barrett formula, which describes the 
saturation of e'(T) at low temperatures due to quantum 
fluctuations. Based on the plasma frequencies of the po- 
lar phonons, the TOl mode was assigned to Slater mode, 
i.e. Ti cations vibration against the oxygen octahcdra, 
and the T02 phonon (near 150 cm -1 ) was assigned to 
Last mode, i.e. vibration of Eu cations against the TiOe 
octahedra.— Temperature dependence of the IR spectra 
indicated coupling (mixing of their eigenvectors) of both 
the modest Moreover, a strong spin coupling with the 
soft mode (SM)is expected microscopically due to a su- 
perexchange between Eu 2+ 4/ spins via the 3d states of 
the Ti 4+ ions^ 

Here we report a study of the epitaxial com- 
pressively strained EuTiC>3 thin films deposited on 
(001) (LaAlO 3 ) . 29 -(SrAl 1/2 Ta 1/2 O 3 ) . 7 i (LSAT) sub- 
strate. Lee et al£ have shown that the maximum ex- 
pected compressive strain of -0.9% in this case is not suf- 
ficient for induction of the ferroelectric and ferromagnetic 
order. Nevertheless, we will show here that the phonon 
properties of EuTiOa/LSAT are strongly influenced: The 
in-plane soft mode frequency is remarkably stiffened in 
comparison to the bulk EuTiC>3 ceramics and moreover 
its IR reflection band is more than five times narrower. It 
enabled us a detailed investigation of the magnetic field 
dependence of this phonon. We will show that shift of 
the lowest-frequency optical phonon with magnetic field 
is mainly responsible for observed magnetodielectric ef- 
fect in EuTiC>3. Moreover, we will report about an anti- 
ferrodistortive (AFD) phase transition appearing in this 
film near 150 K. 



II. EXPERIMENTAL 

Films of two thicknesses (22 and 42 nm) were de- 
posited by reactive molecular-beam epitaxy on the (001) 
LSAT substrate. Details of the deposition were described 
elsewhere J£ The 42 nm film was slightly relaxed, but the 
22 nm film was fully (-0.9%) compressively strained, its 
XRD showing the same in-plane lattice parameter (3.870 
A) as the substrate. The substrates of 10x10x1 mm size 
were provided by R. Uecker and his colleagues from the 
Institute of Crystal Growth, Berlin, Germany. 

The unpolarized IR reflectance spectra were taken us- 
ing a Bruker IFS 113v FTIR spectrometer at temperatu- 
res from 1.8 to 300 K with the resolution of 2 cm" 1 . We 
checked also the polarized spectra at 300 and 10 K, but 
no in-plane anisotropy was observed, confirming that the 
substrate plate as well as the thin films are macroscopi- 
cally optically isotropic in the (001) plane down to 10 K. 
An Optistat CF cryostat (Oxford Instruments) was used 



for cooling the samples without magnetic field down to 
5K. The investigated spectral range (up to 650 cm -1 ) 
was limited by the transparency region of the polyethy- 
lene windows of the cryostat. A helium-cooled Si bolome- 
ter operating at 1.6 K was used as a detector. Custom- 
made superconducting magnetic cryostat was used for IR 
studies at various magnetic fields up to 13 T at tempera- 
tures between 1.8 and 4.2 K. 

Each of the reflectance spectra was evaluated as a two- 
layer optical system^ At first, the bare substrate reflec- 
tivity was measured as a function of temperature and 
carefully fitted using the generalized factorized damped 
harmonic oscillator models 
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where u>TOj arid uiLOj denote the transverse and longi- 
tudinal frequency of the j-th polar phonon, respectively, 
and jTOj and 7 LOj denote their corresponding damping 
constants. s*(uj) is related to the reflectivity R(w) of the 
bulk substrate by 



R(u>) = 



(2) 



The high-frequency permittivity = 5.88 resulting 
from the electronic absorption processes was obtained 
from the room-temperature frequency-independent re- 
flectivity tails above the phonon frequencies and was as- 
sumed to be temperature independent. 

When analyzing the reflectance of the substrate to- 
gether with the film, we kept the earlier fitted bare sub- 
strate parameters fixed at each temperature and only the 
dielectric function of the film was adjusted. For this pur- 
pose, we preferentially used a classical three-parameter 
damped oscillator modet^ 
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since it uses less fitting parameters. This approach is 
well justified, because the damping of LO phonons of the 
film do not influence appreciably the reflectance spec- 
tra. Only the lineshape of the lowest frequency doublet 
was not well reproduced and therefore a coupled damped 
harmonic model was applied at low temperatures (see 
below.) 

The IR reflectance measurements were performed on 
samples with a 1 mm thick substrates which was thick 
enough to avoid the parasitic IR signal from the multiple 
reflection in the substrate. The same samples were af- 
terwards polished to 248±l/im to reach the optimal MW 
resonance near 15 GHz where the samples were measured 
as the composite TE015 dielectric resonators in the shield- 
ing cavity^ using their dielectric resonance. Precisely 
the same thicknesses of bare substrates and substrates 
with the films are needed for an accurate evaluation of 



3 



the in-plane complex permittivity of the films. Each sam- 
ple was measured twice, in the Sigma System M18 tem- 
perature chamber (100-370 K) and in the He-cooled Janis 
close-cycle cryostat (10-370 K). 

High resolution X-ray diffraction (XRD) measure- 
ments were performed on a two-axis diffractometer in 
a Bragg-Brentano geometry (focalization circle with a 
diameter of 50cm). Cu-Ka (wavelength A=1.54056A) 
radiation was emitted from a 18 kW rotating anode gen- 
erator. The out of plane lattice constant of the film as 
well as the substrate were determined with accuracy of 
±0.0003 A. Cryostat was operating between 10 K and 
300 K with an accuracy better than 0.5 K. 

Magnetic properties were measured using a SQUID 
magnetometer (Quantum Design) from 2 to 30 K. Sam- 
ple with the 42 nm thin film was cut to four pieces of 
3.5x4.5x0.248 mm 3 size and all four stacked samples were 
measured simultaneously. 

III. THEORETICAL CONSIDERATIONS 

In parallel with the experimental studies, we also inves- 
tigated the EuTiC>3 films using a first-principles density- 
functional calculations within the spin-polarized gradi- 
ent approximation (GGA). In particular, we have stud- 
ied the effect of the strain and magnetic field on phonon 
frequencies in the bulk and thin films. Our initial guess 
for the low temperature structure of EuTiOa is analo- 
gous to that of the bulk SrTi03 (14 /mem space group). 



In the course of this work, this structure was confirmed 
experimentally^ for the bulk EuTiC>3 and according to 
our calculations, it is also a plausible structure for the 
thin films under the compressive strains. 

For the electronic structure calculations and structural 
relaxations we used projector augmented- wave potentials 
as implemented in the Vienna Ab initio Simulation Pack- 
age (VASP). To account for the strong electron corre- 
lation effects on the /-shell of Eu atoms, we used the 
DFT+U scheme in the Dudarev's approach with an on- 
site Hubbard U and Hund's exchange Jh- For the bulk 
EuTiC>3 of space group Pm3m, U=5.7 eV is used. In 
order to correct for the effect of oxygen octahedron ro- 
tations on the magnetism, the Hubbard parameter is 
increased to U=6.2 eV for EuTiC>3 in the space group 
14/mcm^- Hund's exchange is kept fixed at Jjj=1.0 eV 
for all calculations. Details of the calculations can be 
found in Refs££. 

For the sake of comparison with the experiment, we 
have calculated frequencies of the polar phonon modes in 
the bulk cubic and tetragonal EuTiC>3 as well as in tetrag- 
onal lattices with equatorial lattice parameters fixed to 
those of the SrTi03 and LSAT substrates. In order to 
estimate the frequency changes under strong bias mag- 
netic field, the calculations were carried out both for the 
AFM order (anticipated ground state) and for ferromag- 
netic (FM) order (i.e. for the expected induced magnetic 
order under magnetic field). Resulting mode frequencies 
are listed in Table fl] 



Bulk EuTiOa 


Bulk EuTi0 3 


EuTi0 3 /SrTi0 3 


EuTiOs/LSAT 


EuTiOs/LSAT 


Pm3m 


1 4 /mem 


1 4/ mem 
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Experiment at 1.9 K 
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FM AFM 


FM AFM 


FM AFM 


FM AFM 


FM AFM 


FM AFM 


FM AFM 
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70 77 


125 133 


107 113 


132 138 


77 86 


113 124 


121 125 


104.0, ? 


105.5, 115 


162 165 


158 159 


156 158 


159 161 


147 148 


160 161 


164 166 


? 


162 






251 251 




246 246 




255 255 










422 422 




414 414 




428 428 






549 550 


533 534 


526 527 


536 537 


519 519 


531 532 


532 532 


? 


551 



TABLE I: Polar phonon frequencies in cm -1 of EuTiOa from 
the first principles for the bulk and for films on two different 
substrates, SrTiOs (i.e. no strain) and LSAT (compressive 
strain -0.9%). Experimental results are given in the last col- 
umn for comparison. Only the lowest-frequency phonon was 
studied as dependent on magnetic field. The calculations are 
valid at OK. If the oxygen octahedron rotations are not taken 
into account, the structure of the bulk EuTiOa is Pm3m. If 
the oxygen octahedron rotations are taken into account, the 
space group is (14/mcm), the same both in bulk and thin 
films. The in-plane polarized polar modes have the E u sym- 
metry, whereas the out-of-plane modes have the A 2u symme- 
try. 
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FIG. 1: (Color online) Temperature dependence of the mag- 
netization of 42 nm thick film at different external magnetic 
fields. 



IV. RESULTS 
A. Magnetic properties 

Magnetization curves taken at various external mag- 
netic fields (see Fig. [1]) reveal Tjv=6.3K. It means that 
the -0.9% compressive strain in the film enhances T^ by 
1.0 K because the coupling of spins localized at the 4/ 
levels of Eu cations increases with the reducing lattice 
parameter in the compressively strained films. Note that 
the +1.1% tensile strain in EuTiOa/DyScOa induces the 
ferromagnetic order below 4.2 K<£ On the other hand, it 
was reported that the as-deposited thin EuTi02.86 films 
on SrTiOs (this substrate has exactly the same lattice 
constant 3.905 A as the bulk EuTiOa) exhibit an elon- 
gated c-axis and a ferromagnetic order below 5 K; 21 i 22 
while the same films post annealed at 1000 °C in reducing 
atmosphere exhibit AFM order and relaxed out-of-planc 
lattice strain. 22 



B. Temperature dependence of the phonon spectra 

IR reflectance of the EuTiOa 42 nm film plotted at 
various temperatures is shown in Fig. [3J Comparison 
of the spectra with those of bare LSAT substrate al- 
lows us to distinguish the three TO phonons of the film 
(marked by the arrows). The highest frequency phonon 
is marked T04 (the T03 phonon is silent). The two 
lowest-frequency TO phonons show up softening on cool- 
ing, which is better seen in the inset of Fig. [2j The three 
modes correspond well the three polar modes of the cubic 
EuTiOa, as shown in the Table I. Moreover, a new weak 
mode appears near 112 cm -1 below 100 K. Since the 
42 nm film was not fully strained, we have measured also 
the IR reflectance spectra from a fully (and uniformly) 
strained 22 nm film. Due to the smaller film thickness, 



FIG. 2: (Color online) Temperature dependence of the IR re- 
flectance of the EuTiOs 42 nm film deposited on the LSAT 
substrate. Room-temperature reflectivity of pure LSAT sub- 
strate is also plotted for comparison (dashed-dotted line). 
Phonons from the thin film are marked. Inset shows the en- 
larged low-frequency part of the reflectance spectra, where 
the softening of TOl and T02 phonons is clearly seen. 
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FIG. 3: (Color online) Temperature dependence of the low- 
frequency IR reflectance of fully strained EuTiOs 22 nm film 
deposited on a LSAT substrate. Remarkable TOl phonon 
softening on cooling and activation of a new mode below 150 K 
is clearly seen. 



the intensity of the TOl reflection band is lower (see 
Fig. [3]), but one can clearly see the mode splitting al- 
ready below 150 K. Note that in this case the new mode 
is activated below the TOl phonon frequency, while in 
the thicker film the new mode has higher frequency than 
the TOl mode. Origin of the new mode will be discussed 
in section V. 

Shape of the IR reflectance bands in Fig.[3]is rather un- 
usual and the splitting of the TOl mode cannot be fitted 
well just with a sum of independent oscillators (Eq. [3]). 
Therefore, the coupled oscillator formula was used for the 
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fit of the split modes (their bare parameters are marked ted with Eq. [3l The coupled oscillator formula has the 
with subscripts a and b) and the rest of modes were fit- form22 



Aei a u>% ola (uj% olb - oj 2 + iujjTOib) + Ae lb u>^ olb (uj^ ola - co 2 + iujTOia) - 2iujT^J ' Aei a ^ ola Aei b ^oi b 
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FIG. 5: (Color online) Complex dielectric function spectra of 
22 nm EuTiOs film evaluated at various temperatures from 
the reflectance spectra in Fig. [3] The inset shows the temper- 
ature dependence of the imaginary coupling constant between 
two lowest-frequency modes. 



FIG. 4: (Color online) a) Temperature dependence of the 
phonon frequencies for the EuTiOs films of two thicknesses 
42 and 22 nm below 170 cm" 1 . Stiffening of the TOI and T02 
phonon frequencies in the more strained 22 nm film is seen. 
The new mode in the 22 nm film activates by 50 K higher than 
in the thicker film. The inset shows the temperature depen- 
dence of the new mode frequency in expanded scale down to 
1.9 K. Mode hardening below the Neel temperature is clearly 
seen, b) Temperature dependence of the oscillator strengths 
of all the modes seen below 200 cm -1 . One can see that the 
new TOla mode gains strength from the TOlb mode. 



using the imaginary coupling constant T model. The rest 
of symbols have the same meaning as in Eq. [3] 

Fig. [T^i shows the temperature evolution of the po- 
lar phonon frequencies below 170 cm -1 obtained from 
the reflectance fits of both films. At 10 K the TO 



phonon frequencies are 83, 151 and 533 cm -1 for the 
bulk ceramics^ while our 22 nm film shows phonons at 
115, 162, 551 cm -1 and the additional new mode at 105 
cm -1 . There is a strong phonon frequency stiffening in 
the films due to the compressive strain. Note also that 
the phonon frequencies are higher in the 22 nm film than 
in the 42 nm one (see Fig. |4^) , obviously due to the higher 
strain in the former case. 

Complex permittivity spectra evaluated from the IR 
reflectance are shown in Fig. [5j Stiffening of all in-plane 
polar phonons is responsible for a lower in-plane permit- 
tivity e' in the films compared with the bulk sample. 
At 10 K the static permittivity reaches the value of 146 
compared to 405 in the single crystal.— Here it should 
be stressed that compressive strain reduces the in-plane 
permittivity, but the out-of-plane permittivity should be 
enhanced, because the TOI phonon component polarized 
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perpendicularly to the film plane should be softer. Unfor- 
tunately, we cannot see this phonon in our near-normal 
reflectivity geometry. 

In our fit we used Aei a =0, i.e. we suppose that the 
bare lowest-frequency mode has zero dielectric strength. 
The imaginary coupling does not shift the phonon fre- 
quencies, it merely deforms the spectral line shape and 
transfers the oscillator strength Aeuitoj fr° m TO lb to 
the new TOla mode (see Fig. HJd). Dramatic increase of 
the coupling constant T is seen on cooling below 150 K 
in inset of Fig. [5] This causes the increase of its strength 
on cooling. Also the temperature dependences of the 
TOlb and T02 mode oscillator strengths are remark- 
able. They give evidence on mixing of eigenvectors of 
both the modes. 




50 100 150 200 250 300 350 400 
Temperature (K) 



C. Microwave dielectric measurements 



The non-monotonous dependence of the oscillator 
strengths of TOlb and T02 modes (see Fig.|3b) could be 
related to an AFD phase transition involving anti-phase 
tilts of the oxygen octahedra (a°a°c~ in Glazer nota- 
tion). This phase transition occurs in the bulk EuTi03 
near room temperatur e 8 ' 10 , but in the strained films it 
can be expected at lower temperatures (e.g. near 180 K 
in the 1% tensile strained EuTi0 3 /DySc0 3 ) 3 . Also the 
TOla -TOlb coupling could be related to this transition, 
if the distorted phase would be ferroelectric. Therefore, 
we have investigated the temperature dependence of the 
in-plane MW permittivity of the film and of the sub- 
strate itself. Complex permittivity of the film measured 
near 15 GHz shows a dielectric anomaly near 120 K (see 
Fig. [5]) reminding a ferroelectric phase transition. Nev- 
ertheless, the second harmonic generation (SHG) mea- 
surements of the EuTi0 3 /LSAT film did not reveal any 
acentric phase down to 2 K. 3 Moreover, from theory we 
know that the spontaneous polarization (and the main 
dielectric anomaly) should occur in the [001] direction 
in the compressively strained films,— whereas we mea- 
sure only the in-plane dielectric response. It is suspicious 
that the low temperature MW g' is higher than the sum 
of phonon contributions to e' (see Fig. HJ) and that the 
maximum of e' occurs 30 K below the temperature where 
the new mode activates in the IR spectra. Since we do 
not see any phonon anomaly in the IR spectra explain- 
ing the dielectric anomaly near 120 K, it seems that the 
MW dielectric anomaly has its origin in some dielectric 
dispersion below the phonon frequencies, most probably 
due to diffusion of charged defects (e.g. oxygen vacan- 
cies). Such defects play the main role in the conductivity 
of the films, which deteriorates the dielectric studies of 
EuTi0 3 /LSAT in the kHz region. 



FIG. 6: (Color online) Temperature dependence of the com- 
plex permittivity in the 22 nm EuTiC>3 film (left scale) and 
the LSAT substrate (right scale) measured at 15.75 GHz. The 
thin film MW data were obtained using two low-temperature 
devices (LN2 chamber and He refrigerator) . Static permittiv- 
ity calculated from the phonon contributions (solid dots) is 
shown for comparison. 



D. Magnetodielectric effect and tuning of the 
phonon frequency by magnetic field 

Bulk EuTi0 3 has a high static s' ~ 400 and low SM 
frequency. Therefore it shows a high (7%) tunability of e' 
by the magnetic field 5,— and an unusually high 7 cm -1 
change of the SM TOl frequency with B is expected 
theoretically— Such a shift was not confirmed in our IR 
spectra of the bulk ceramics, presumable because the 
TOl reflection band is too broad (almost 100 cm -1 )M- 
However, the sharper TOl reflectance band in compres- 
sively strained EuTi0 3 films (see Figs. [2] and |3j) appears 
to be more promising for such observation, although the 
expected frequency shift (see Table [TJ is somewhat lower 
(ps 4cm~ 1 ). 

Actually, we have observed an almost 2 cm -1 decrease 
of the lowest-frequency phonon with B up to 10 Tesla 
(see Figs. [7] and [5]). The phonon frequency shifts could 
be determined with a high accuracy of 0.1 cm -1 thanks 
to the sharp feature of the reflectance band. The absolute 
changes of the phonon frequencies on B were similar and 
reproducible in both films, even though the origin of the 
sharp phonon seen in Figs. 2 and 3 appears to be differ- 
ent. It corresponds predominantly to TOl and an optical 
phonon from the Brillouin zone edge (see the discussion 
in the next section) in the 42 and 22 nm film, respectively. 
The higher frequency component of the TOl reflectance 
band is weaker and broader, therefore it is not possible 
to evaluate its change with B. Magnetodielectric effect, 
i.e. the relative change of the static permittivity with 
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FIG. 7: (Color online) IR reflectance spectra of the 22 nm 
thin film at 1.9 K taken at various magnetic fields. The shift 
of phonon frequency is clearly seen. 



magnetic field B (Ae' (B) / e' (0)) obtained from the fits 
of reflectance spectra is plotted in Fig. |9] As expected, 
one can see that the change of Ae'(B)/e'(0) in the thin 
film is almost three times smaller than that in the single 
crystal,— because the phonons are stiffened in the films. 
The change of Ae'(B) /e'(0) appears to be larger at 1.9 K 
than at 4.2 K (Fig. [HJ) - It is caused by the larger phonon 
frequency shift with B at 1.9 K (see Fig. [8]), since also 
the shift of the phonon frequency on cooling below Tjv 
is larger at 1.9 K than at 4.2 K (see inset of Fig. |4j. It 
corresponds also to larger relative change of the permit- 
tivity at 2K than at 4.2 K reported in single crystal^ 
and ceramics.— Here it should be noted that the mag- 
netic field influence on the A^ u SM component polar- 
ized perpendicular to the film plane (unattainable in our 
experiment) should be higher than the influence on the 
phonons polarized in the film plane. This is shown also 
in Table I. 

Phonon frequencies obtained from the first principles 
calculation in the AFM (i.e. without applied B) and 
in the ferromagnetic phase (i.e. in a high B) are sum- 
marized in Table HI One can see that mainly the lowest 
frequency phonons shift down with the magnetic field. 
Theoretical frequency shift is twice larger in Table fl] than 
in our experiment (Fig. [3]), because the calculations were 
performed at conditions of K, while the experiment was 
done at «2K and the shift should increase on cooling 
below Tn- In this way we confirmed that the magne- 
todielectric effect in EuTiC>3 is due to the change of the 
polar phonon frequency with magnetic field. 

It is well known from the literature that it is rather easy 
to change the phonon frequencies by hydrostatic pressure 
or external electric field. This effect is most remarkable in 
ferroelectrics, in which the ferroelectric SM is highly sen- 
sitive to boundary conditions. Electric field dependence 
of the ferroelectric SM was first-time observed in Raman 
scattering of SrTi03 ) 25 i 26 and the hydrostatic pressure 
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ative changes of static permittivity obtained from the fits of 
IR reflectance of the 22 nm thin film. 



dependences of the SM were reviewed by Samara and 
Peercy<21 

On the other hand, direct observations of phonon fre- 
quency changes on application of external magnetic field 
B are very scarce in literature and most of magnetic ma- 
terials do not exhibit measurable effects. Ruf et ali 2 ^ 
observed 1.5 cm" 1 change of one phonon line with B 
in the Raman scattering spectra of YBa2Cus07. This 
phonon frequency was sensitive to the superconductive 
phase transition and its changes were observed slightly 
below T c under the field of 6-12 T. Sushkov et al^ re- 
ported on small B-dependent line-shape changes in the 
IR absorption spectra of molecular magnet Mni2-acetate, 
but a very high magnetic field of 30 T was necessary to get 
a measurable effect. Multiferroics exhibit relatively high 
tunability of e' so that also larger changes of phonon fre- 
quencies with B could be expected. This was confirmed 
in GdMn0 3 and DyMn 2 5 multiferroicsi^ML In the for- 
mer case, a phonon frequency change by less than 1 cm" 1 
was observed due to a phonon coupling with magnon^ 
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while in the latter case a high field of 18 T was neces- 
sary to see a ^2 cm -1 change of phonon frequency^ In 
our EuTiOa film the phonon frequency change due to 
B is comparable to the previous case, but a much lower 
magnetic field was needed. 

It should be stressed that the magnetodielectric effect 
in EuTiOs is not caused by a linear magnetoelectric cou- 
pling, since bilinear terms aijHiEj (where Hi and £j 
are components of the magnetic and electric field, respec- 
tively) are not allowed in the thermodynamic potential. - 
Nevertheless, higher order magnetoelectric effects, ac- 
counted for by the PijkHiHjE^ and Sij^iHiHjE^Ei terms 
in the thermodynamic potential, are permitted. The 
strong quadratic and bi-quadratic magnetoelectric effects 
were actually experimentally confirmed in bulk EuTi03 
by Shvartsman et alJ^ 

Dzyaloshinskii recently proposed a splitting of polar 
phonons in ferroelectrics due to an external magnetic 
field.' 52 We do not see any splitting of phonons in the 
magnetic field, but we observe a 50% increase of the SM 
damping with magnetic field (see the broadening of re- 
flection band with magnetic field in Fig. [7]). At 1.9 K its 
damping is 4.2 cm -1 at 5=0 T, while at 5=13 T it is 
6 cm -1 . Broadening of the phonon linewidth could be 
a consequence of a small phonon splitting, which is not 
resolved. 



V. DISCUSSION 

Let us now try to explain the origin of activation of 
the new mode near the TOl phonon frequency in the IR 
spectra of EuTi0 3 /LSAT (see Fig. gj. We consider two 
possibilities: 1) the doubly degenerate E u component of 
the TOl phonon probed in our experiment splits due 
to anisotropic in-plane strain, which is created due to 
a phase transition in the LSAT substrate. 2) the new 
mode originates from the Brillouin zone (BZ) edge due 
to an AFD phase transition (with multiplication of the 
primitive unit cell and BZ folding). 

Let us discuss the first possibility. The in-plane polar- 
ized modes can split if the tetragonal symmetry of film re- 
duces to orthorhombic or lower- symmetry structure. It is 
possible only in the case, when the in-plane strain would 
be anisotropic, i.e. the LSAT substrate would not be 
cubic. Indeed, it has been reported that the LSAT sub- 
strate exhibits a small distortion from cubic symmetry 
near 150 K and its low-temperature phase has tetragonal 
or orthorhombic symmetry^ We have measured temper- 
ature dependence of the LSAT lattice constant and it re- 
ally exhibits a small hump near 150 K (see Fig. [T0|) and 
the same anomaly is revealed also in the out-of plane 
EuTi0 3 film lattice parameter. Moreover, the struc- 
tural phase transition in LSAT can be recognized also 
in the microwave dielectric permittivity, where a small 
anomaly was observed near 150 K (see Fig. Never- 
theless, in case of substrate-induced anisotropic strain 
the split modes should have different symmetries and 
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FIG. 10: (Color online) Temperature dependence of the out-of 
plane lattice parameter of the EuTiC>3 film (thickness 22 nm) 
and of the LSAT substrate. 



therefore they cannot couple, as observed in our spectra. 
Moreover, the possible anisotropic strain from the non- 
cubic LSAT substrate should as well split the phonons 
in SrTi0 3 /LSAT (SrTi0 3 and EuTi0 3 have the same 
lattice constant and therefore the strain in both films 
is the same), but this was not detected^ On the other 
hand, we have revealed the same new phonon activated 
in IR spectra near 150 K in EuTi0 3 films deposited on 
NdGa0 3 and LaA10 3 . Compressive strains in both films 
are slightly relaxed and have the value of about 0.9%, 10 , 
i.e. close to EuTi0 3 /LSAT. Since the NdGa0 3 and 
LaA10 3 substrates do not exhibit any structural phase 
transition, the new mode observed in all three films can- 
not be activated due to anisotropic strain. 

Let us discuss the second possibility that the new mode 
originates from the BZ edge. It can be expected that the 
AFD phase transition evidenced in the bulk EuTiO^ - — 
and tensile strained thin film EuTi0 3 /DySc0 3 (see sup- 
plement of Refi^) should occur in the compressive thin 
films as well. The tetragonal I A/ mem symmetry reported 
for the bulk can be also well compatible with the thin 
film epitaxial geometry. Temperature of the phase tran- 
sition may be considerably shifted in comparison to the 
bulk, of course. In the AFD phase two new E u phonons 
near 255 and 430 cm -1 stemming from the R-point of 
the BZ edge become active in the IR spectra, as it can 
be seen from the Table I. However, the frequencies are 
much higher than that of our new mode seen between 
110 and 105 cm -1 . 

It is known that two structural soft modes (SMs) from 
the R-point of the BZ edge can be activated in low- 
frequency Raman spectra below the AFD phase tran- 
sition in SrTi0 3 (T C =105K),2I The AFD structure of 
SrTi0 3 is the same 1 4/ mem as in EuTi0 3 and the two 
AFD SMs should be IR inactive in both materials, if 
their structure remains centrosymmetric. However, these 
AFD SMs were also discovered in the low-temperature 
IR spectra of bulk SrTi0 3 ceramics^ and polycrystalline 
films' Their activation was explained by polar distor- 



tion near the grain boundaries j 38 i 39 Activation of 
AFD SM was also observed in tensile strained SrT 
and EUT1O3 films on DyScC>3 substrates, since the 
sile strain induces a ferroelectricity in both films,— 
Our new mode seen near 105 cm -1 does not remind 
AFD SM, because it does not harden below the assm 
AFD transition temperature near 150 K. Nevetheles; 
can be an optical mode from the BZ edge. Theoi 
cal calculations^ show that the lowest-frequency opt 
branch exhibits almost no q dispersion in BZ, there 
TOl mode (q=0) and the optical mode from the BZ c 
have almost the same frequency. AFD SMs have pi 
ably much lower frequencies, therefore they are not 
solved in our spectra. Nevertheless, both AFD SMs 
optical mode from BZ boundary can be activated a 
folding of BZ zone only if the crystal structure beco 
at least locally broken. Even if we know that the los 
the inversion center was not evidenced by the SHG 
nal in EuTiOa/LSAT, 3 , it is known that the IR probe is 
sensitive on local breaking of the symmetry known e.g. 
from polar nanoclusters in relaxor ferroelectrics42 

Here we propose a new possible mechanism of po- 
lar nanocluster appearance in the EuTiOs as the con- 
sequence of AFD phase transition and oxygen vacancies, 
which are always present in the lattice. Let us imag- 
ine that we have oxygen vacancies in the Ti02 plane 
of the perovskite structure (see Fig. [TT]). Then, in the 
ideal perovskite structure without octahedra tilting the 
local environment around the vacancy is antipolar: two 
neighboring Ti atoms are equally displaced in opposite 
directions. In a structure which allows oxygen octahedra 
tilting, the displacements of Ti atoms are antipolar in one 
direction, but with uncompensated dipoles in orthogonal 
direction (see Fig. [Til) . In such case the local inversion 
symmetry is broken. We have calculated the structure 
in a 40-atomic simulation cell (i.e. EuTi02.875 instead of 
EuTiOa) around this defect and found that it has caused 
Ti displacements not only around the vacancy, but also 
in other Ti's in the vacancy neighborhood. Of course, 
the displacements of the neighboring Ti are the largest. 
Also, we found that the formation energy of the O va- 
cancy in the Ti02 plane is by 0.03 eV lower than of the 
vacancy in the EuO plane (i.e. apical oxygen). 

The IR inactive folded modes can become IR active 
also if we assume AFM short-range order even if it re- 
mains centrosymmetric. Our density-functional calcula- 
tions for an auxiliary /42m structure compatible with 
AFM order yielded a new IR active phonon (originally 
optical one from the BZ edge) with frequency that is 
12 cm -1 lower than the TOl mode. The strength of this 
mode strongly depends on the conditions under which 
it is calculated. We have used finite displacements ap- 
proach to calculate the force constant. For 0.025 A dis- 
placement from their equilibrium position, the mode- 
plasma frequency V Aeujto of the new mode is about 
270 cm -1 (for comparison the TOl mode has the plasma 
frequency around 1040 cm -1 ). If the displacement is only 
0.0027 A, then the calculated mode-plasma frequency is 
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FIG. 11: (Color online) Schematic plot of the perovskite crys- 
tal structure in the view along [001] axis in the a) untilted and 
b) tilted AFD phase. Squares are marks for the oxygen octa- 
hedra, Eu cations are not shown. Orientation of the electric 
dipole moment created by the Ti displacement in the vicinity 
of oxygen vacancy in the tilted phase is marked by a large red 
arrow. 

only 24 cm" 1 . This is a sort of probe of anharmonicity 
and it points to a rather strong anharmonicity of this 
particular mode. This mechanism of the mode activa- 
tion requires a short-range magnetic order at least up to 
150 K. Recent not yet published /iSR experiment^ 3 - per- 
formed on bulk EuTiOa revealed short-range magnetic 
correlations up to 300 K, which is rather surprising, be- 
cause the AFM critical temperature is only 5.3 K. 

Finally we can conclude that the new mode near 
105 cm -1 is of the same symmetry as the TOl mode. 
It probably stems from lowest-frequency optical phonon 
branch and it is activated from the BZ edge due to the 
AFD phase transition. Its activation in IR spectra re- 
quires local reduction of the symmetry, which can occur 
due to creation of polar or magnetic nanoclusters. 

VI. CONCLUSION 

Due to the compressive strain, all polar phonons po- 
larized in the plane of the EuTi03 films are consider- 
ably stiffened in comparison to the bulk ceramics. Below 
150 K a new IR active mode was revealed near 110 cm -1 . 
We understand it as optical phonon from the BZ edge, 
which activates in the IR spectra below the AFD phase 
transition under circumstance of presence of polar and/or 
magnetic nanoclusters. Magnetization data show that 
the compressive strain enhances the Neel AFM tem- 
perature by 1 K due to the strain enhanced superex- 
change interaction. IR reflectance studies in magnetic 
field revealed a striking tuning of the lowest phonon fre- 
quency with the magnetic field. This is responsible for 
the magnetodiclcctric effect reported previously in the 
bulk EuTi03. n Observed change of the polar phonon 
frequency with the magnetic field is in agreement with 
our theoretical calculations. Similar IR experiments with 
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magnetic field are rather rare in the literature, because 
the phonon frequency changes are usually under spec- 
tral resolution in materials without a strong spin-phonon 
coupling. Our study shows that the IR reflectance mea- 
surements of a rather thin film can be more sensitive to 
the magnetic field (in the case of suitable non-conducting 
substrate) than the reflectivity spectra of the bulk sam- 
ples or transmission spectra of the thin films. So, the 
IR reflectance spectroscopy is a very promising tool for a 
study of the magnetodielectric effect as well as of struc- 
tural phase transitions in thin dielectric films deposited 
on electrode-less nonconducting substrates. 
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